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Abstract 

Quantum field theory is applied to study the interaction of an electron 
plasma with an intense neutrino flux. A connection is established between the 
field theory results and classical kinetic theory. The dispersion relation and 
damping rate of the plasma longitudinal waves are derived in the presence of 
neutrinos. It is shown that Supernova neutrinos are never collimated enough 
to cause non-linear effects associated with a neutrino resonance. They only 
induce neutrino Landau damping, linearly proportional to the neutrino flux 
and Gp. 
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1 Introduction 



Bingham et al. |l| have studied the interaction of a neutrino beam with an electron 
plasma and concluded that the neutrino fluxes produced in Supernovae are intense 
enough to cause plasma instabilities with large growth rates. If true this would 
provide a physical mechanism of energy transfer from the neutrinos to the plasma 
that might explain the Supernova explosions. The interaction between neutrinos 
and plasma was described with a ponderomotive force acting on the electrons and 
a neutrino wave function obeying a naive Klein-Gordon equation modified with 
a matter induced external potential. This non-standard treatment was not well 
established from the Standard Model of electroweak interactions and originated 
some controversy P, |!J. More recently classic kinetic theory was applied to 
study the neutrino-plasma system where both neutrino and electron particles suffer 
each own ponderomotive force. Again the lepton spin and chiral structure of the 
weak interactions remain unnoticed. The dispersion relation derived for the plasma 
waves differed from the one previously obtained in |l|. Yet, the authors reiterated 
the claim that the neutrinos produce non-linear effects for certain resonant modes 
of plasma waves causing instabilities with large growth rates proportional, not to 
Gpn u , but to a smaller power of this quantity. 

The aim of the present work presented here is to obtain a formulation of the 
problem based on field theory and Standard Model of electroweak interactions, an 
effort initiated in ||. We confine to isotropic plasmas and longitudinal photon 
excitations (also called plasmons or Langmuir waves). The Cerenkov emission of 
longitudinal photons by massless Standard Model neutrinos in an isotropic plasma 
has been studied along the years || [7|, ||, |9], |ID| . This is a single neutrino decay 
but not a collective neutrino process. Hardy and Melrose |J gave one step more 
by extending the work to spontaneous and stimulated emission and absorption of 
plasmons to study the so-called kinetic instabilities (see also [[11]). They obtained 



an expression for the decay (growth) rate of the plasma waves induced by a neutrino 
flux. However, because the derivation was based on single neutrino processes the 
result is necessarily proportional to G| excluding apriori any possible non-linear 
effects. In the present paper we derive the neutrino contribution to the photon self- 
energy and obtain a modified dispersion relation for the longitudinal waves. This 
allows one to study not only kinetic instabilities but also the possible existence of 
hydrodynamic instabilities and non-linear phenomena. 

Next section we calculate the neutrino contribution to the electromagnetic po- 
larization tensor and the dispersion relation of longitudinal photons. In section 3 we 
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establish a relationship with kinetic theory and in section 4 the dispersion relation 
and nature of neutrino induced instabilities are analyzed in detail. In the last section 
we summarize the main results. 



2 Electromagnetic polarization and longitudinal 

waves 

In a medium the Maxwell equations are modified by polarization effects. The elec- 
tromagnetic waves obey the following equation: 



7T 







where e is the wave polarization vector, k the linear momentum and is the 
polarization tensor. The purely electromagnetic contribution, vr^, is well known 
for an electron plasma in first order of approximation and can be identified with 
the Feynman diagram of Fig. [I]. At finite temperature the one-particle propaga- 
tors possess additional terms related to the one-particle distribution functions of 
the background matter. In the real-time formalism [13], [TJJ], adopted here, the 
electron propagator in a homogeneous and unpolarized plasma is given by 



(p 1 + m) 



P 



2n5(p 2 - m 2 ) (9(p°)f e (p) + 9(-p°)f g (-p) 



(2) 



where f e , f s are the electron and positron distribution functions respectively. 

Gauge invariance implies that e is defined up to a vector proportional to k 
and kf/ir^ = = n^ v k v . On the other hand, for an isotropic, homogeneous and 
unpolarized plasma the tensor is symmetric and can be written in terms of the 
metric tensor, momentum k and vector u defining the time direction in the plasma 
rest frame fl3 |. As a result, 

-tt^m (3) 
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Figure 1: Electron plasma contribution to the photon self-energy, m'- 
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Figure 2: Neutrino-plasma contribution to the photon self-energy, iir^ . 

in the plasma frame (u = k°) and the wave Eqs. (|l|) admit a purely electrostatic 
solution, = (1, 0) = u M in the Coulomb gauge. The dispersion relation of these 
longitudinal waves is 

P + vr° M = 0. (4) 

In the low energy limit of the Standard Model of electroweak interactions the 
electron and neutrino interactions (both charged and neutral currents) are described 
by the effective Lagrangian 

Ant = eA fl e^e - \Pl G F (ulI^l) e^(c' v - c' A 7 5 )e , (5) 

where e denotes the positron charge (a = e 2 /47r), Gf the Fermi constant and 

c' v = c' A + 2 sin 2 6 W , (6) 

with c' A equal to +1/2 for v e and —1/2 for u^, v T . The neutrino flux contributes 
to the electromagnetic polarization through the diagram of Fig. The diagrams 
with fermion self-energy corrections either in vacuum or in matter will be neglected 
as well as the weak interactions between electrons or nucleons of the medium. The 
expression of the Vl propagators is similar to the electron propagator of Eq. (fj): 



G v {p) = 



- 2 - 2vr<V) (0(pP)/„(p) + 0(-p°)/„(-p 



(7) 



The only differences are that the v masses are taken to be zero, there is only one 
spin degree of freedom, vl or Ur, and the neutrino and anti-neutrino distribution 
functions, f v and f p , are not thermal, as they move along a privileged direction. 
The neutrino loop gives a tensor, 

= - / ^-jT{ la G u (p + k)^G u (p)} , (8) 

whose vacuum contribution is to be ignored. Only the contributions from the elec- 
tron and neutrino loops that are linear in the respective particle densities will be 
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Figure 3: Neutrino electromagnetic coupling induced by an electron plasma. 

retained. The electron loops are directly related to the neutrino electromagnetic 
coupling. Writing the vvj vertex in Fig. |^ as 

-<7*^r""(A0, (9) 
where k is the incoming photon momentum, the diagram of Fig. [2] is given by 

i7^ = ir^(-fc)7^(fc)i^(fc). (io) 

The ui/j vertex separates in a pseudo-tensor proportional to c' A and a tensor that is 
proportional to TTwyr, as follows M, K 



T^(k) = -V2G F (c'ynZi ~ c> 5 e^k a up) . (11 



For future reference let us write the expressions of the susceptibility tensors, 

ixv o 2 f d3 P fe + fe (J^Y + ~ k-pg^k-p - Ifp^f 

*EM- & J (27r)3 ^ (k-pf-(k*/2f ' K ] 



71 » 



d 3 P fu + fp {^p v + p^k v -k-p g^)k-p - k 2 p»p v 



N J (2tt) 3 E v {k-p) 2 - (P/2) 2 

e^k a p p J {2 J )3 fu Eu fu {k . p) 2_ {k 2 /2) 2 ■ ( 13 ) 

If the neutrino distributions were isotropic the polarization tensor tt^ would 
satisfy a relation like (Q) and the longitudinal photons could still be described with 
a scalar potential i.e., = (1,0) in the plasma rest frame. That is not strictly the 
case in the presence of a neutrino flux but since the plasma is isotropic, as considered 
here, the vector components e l are expected to be suppressed with respect to e° by 
a factor of G F k A . We now prove that the components e l are in fact suppressed and 
their contribution to the dispersion relation only comes at G F order. First, note that 
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the assumed plasma isotropy implies not only the relation (H) but also the structure 



u 2 k l ti 



m 

4 = -KTPij + , (14) 

where P is the projector over the plane orthogonal to k, 

Pij = $ij fT ■ ( 15 ) 

k 2 

Due to gauge invariance only three of the wave Eqs. ([]]) are linearly independent. 
Working in the Coulomb gauge, k-e = 0, the application of the projector P on the 
left of the last three Eqs. (P gives 

(k 2 -n T )e l = P tJ ^e u . (16) 

The factor k^—n? vanishes for transverse photons but not for longitudinal waves: in a 
non-relativistic gas, ttt is equal to u 2 the square plasma frequency, and k 2 = u 2 — k 2 . 
The above equation can be used to calculate e % in a iterative way, as follows: 

e i = Pi j 7^—e° + --- • (17) 
k 2 — n T 

It is clear that e % are suppressed by G| with respect to e°. When substituted in the 
first of the wave Eqs. (]]]), 



it becomes evident that the vector components e l only reflect in the dispersion 
relation at Gp order and can thus be safely neglected. The dispersion relation of 
longitudinal waves is so 

k 2 + C + 74° = (19) 

at Gp order. It remains to evaluate tt^ u and 7r™. 

The electroweak interactions conserve the electric charge and lepton numbers L e , 
Lfj,, L T . This implies that at energies considered here, far below the muon mass, the 
electron, v e , and v T numbers are separately conserved which materializes in the 
conservation laws 

Vr* = = o = <K = r^K . (20) 

Using repeatedly these relations together with Eq. (|3D one obtains from Eqs. ( [10| - |TTD 

^2^/2 „ / , ,2 \ 2 



00 2GpCy / nn \2 ( _ U! 
7T" 



- — — (*) ( 1 - j 2 ) <° • (21) 
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The outcome is that the c' A coupling does not contribute to the dispersion relation of 
longitudinal waves in an isotropic electron plasma. On the other hand, c'y is much 
smaller for u^, v T (c' v ~ —0.04) than for u e (c' v ~ 0.96) so only u e and u e will be 
taken in consideration. 

The expressions of the response functions are 

""em = -2e 2 I ^ (/ e + h) E e J 2 ~ {t E^ , (22) 



(27r) 3 we JSJ e {k . pe) 2_ {k 2 /2 ) 



- - l^-AU + h) K n *J k '%* , (23) 



(27r)3 w ■>»' "(k- Pl/ ) 2 -(k 2 /2) 2 
where E and v denote in each case the energy and velocity of the particle. The 
factor of 2 in front of the electromagnetic susceptibility stands for the number of 
electron spin states. The normalization of the distribution functions is fixed by the 
propagators (Q) and (0) as follows: the number of particles (e~, e + , v or v) per unity 
of volume and spin degree of freedom is 

-/<&'■ ^ 

Next section we seek a connection between field theory and classical kinetic 
theory in the limit of small frequencies and wavelengths. In section 4 we study the 
dispersion relation and possible neutrino induced instabilities. 



3 Kinetic Theory 

The Feynman diagrams in Figs. 1 and 2 are a consequence of a certain field dynamics 
that couples boson fields with fermion densities. The electroweak interactions in 
particular only involve vector fields and currents. The ones relevant here are the 
electromagnetic field and current, and J^ M , plus the neutrino and weak electron 
currents defined as 

J& = v L -fv L , (25) 
J^ e = V2G F eY(c' v -c' Al5 )e. (26) 

Behind those Feynman diagrams there is a set of equations relating the current 
fluctuations, J^(k) in momentum space, to each other: 

^ M = -<u{k)A v + J Na T^(-k), (27) 
J Na = J We , (28) 

J^ e = T^{k)A y + ^ e {k)J Nu . (29) 
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Here, 7Tw e is an electron loop, suppressed by G F , that couples the weak electron and 
neutrino currents. One obtains from this a relation between the electromagnetic 
current and A^. Using matrix notation and the definition r ,fMX (k) = r ati (—k) one 
gets 

Jem = -(7T E M + r / 7r A r(l + 7r We 7r iV )~ 1 r)v4 (30) 
~ -(7T EM + T'tt n T) A , (31) 

where the corrections of G F order associated with 7Tw e are neglected in the last 
equation. When this electromagnetic current is put in the Maxwell equations, 

(_fc2^ + k ^ K = ? (32) 

they give rise to the wave equations ([[]) with a polarization tensor given by the 
Feynman diagrams of Figs. 1 and 2. 

We have seen in the previous section that for an isotropic plasma the c' A coupling 
and electron axial-current do not play a significant role. If that current is dropped 
out of Eqs. (p7)-f29|) they reduce to 

Jem = -^em (~ eA u + V2G F Cy J Nl )j , (33) 

J Na = \i£p (V2G F c' v jP M ) . (34) 

which now have a classical counterpart in the sense that the fermion polarization 
does not appear explicitly. This is the kind of relations also obtained in classical 



kinetic theory, a framework that has been used by some authors |4|, [17]]. That fact 
motivated us to formulate our own classic theory in order to better understand the 
differences between the results based on quantum field theory and their works. 

In kinetic theory a system is described by distribution functions on the phase 
space of the particles, in particular, the single particle distribution functions f(t, x,p). 
In low dense plasmas the collisions are less important than the collective interactions 
and the time evolution of the distribution functions is given by the Vlasov equations 
II, 

?L + v-— + F-— = (35) 
dt dx dp 

The functions velocity, v, and force, F, have to be specified for each of the particles 
species. 

If one ignores the electron polarization and c' A coupling, the effective Lagrangian 
of Eq. d5|) reduces to 

£ int = e A^ - V2G F c' v J Nfl J? . (36) 
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It clearly admits a classical limit where the neutrino vector current is equal to 
the difference between the v and v current densities, J% = j£ —jp, and J£ = — jg 
is the analogous current for electrons and positrons. The next step is to write down 
the interaction Lagrangian for a classical electron or neutrino particle respectively 



L e = (eA' i -V2G F </ v j£)x ll , (37) 
L v = -(V2G F c' v J^ Xfl . (38) 

They are symmetric to the positron and anti-neutrino Lagrangians. Notice that the 
vector current Jj^ couples to an electron particle in exactly the same way as the elec- 
tromagnetic potential and in turn the neutrinos interact with a vector potential 
as well, proportional to Jg. Therefore, the electroweak forces are a straightforward 
generalization of the electromagnetic Lorentz force || . The total force acting on an 
electron is 

F e = -e(E + v e A B) + V2G F c' v (E N + v e A B N ) (39) 
with weak-electric and weak-magnetic fields given by 

d f 

En = - V J° N - , B N = V A J N . (40) 

The positron force is —F e . In a similar fashion, the neutrinos suffer a weak force 

F v = V2G F c' v (E e +v u A B e ) , (41) 
(— F v for anti-neutrinos) where 

df - 

E e = - V J e ° - , B e = V A J e . (42) 

Above, the velocity and linear momentum are related to each other by v = p/ y/p 2 + m 1 
with a zero mass in the neutrino case. 

These weak forces differ from the ones employed before |], [17[]- in particular, 



the ponderomotive forces considered in |l|, Q only contain the terms proportional 
to the gradients, V ' J% and VJ° , of the neutrino and electron densities but not the 
terms that go with the vector currents Jn and J e . In the case of the longitudinal 
waves, the weak-magnetic forces are suppressed by an additional power of G F but 
the time derivatives of Jn and J e still contribute at the same level as the density 
gradients. This fact will manifest in the dispersion relation itself. 
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The proper modes of a system are usually investigated by expanding the distri- 
bution function around a static and uniform function f° i.e., 

f(x,p) = f°(p)+5f(x,p). (43) 



The Vlasov equations (|35|) are then approximated to the linearized form 



' d d\ - df° 



which become [|18 



df° 

-i(uj-k-v)5f(k,p) + F(k,p)~ = (45) 

op 

— * 

after Fourier analysis, where /c M = (u;, k) denote the frequency and wave vector 
components. It is convenient to write this in a relativistic covariant notation. In 
terms of 

iT=P°^, (46) 
that transforms as a 4-vector, the equation above reads as 

df° 

- ik-p 5f(k,p) + F^k,p)^ = . (47) 

It does not matter whether f° depends on p° or it is evaluated on-shell (p 2 = m 2 ) 
because F£ d/dp^ is a total derivative when applied on functions that do not depend 
on the space coordinates. 

Under an external perturbation the current densities 

f(x) = J T0^f( x >P)P" ( 48 ) 
(E is the kinetic energy) suffer fluctuations equal to 

3 [ ] ~ % i (2n)3Edp» hp ' 1 J 

To compare with the field theory results it is convenient to integrate by parts ob- 
taining 
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It must be kept in mind that here dF%/dp v identifies with EOF /dp after setting 
the on-shell condition p° = E. In the cases of interest to us the particles interact 
with a vector potential, V^, and the 4- force is 

F£ = {jPV v - drV)p v . (51) 

Eqs. ( ^9| , |50|) yield a linear relation, = —ti^ v V V) with susceptibility given by 

^ = /^g *y-*™-> (52, 



(2nfEdp a k-p 



or, after integrating by parts, 



r d 3 p f° (k»p v + p^k u - k-pg^k-p- k 2 p»p u 

n [/J ""i(27p¥ Jj^pf • (53) 

For instance, to obtain the electromagnetic polarization tensor it suffices to make 
Jg M = qj^ 1 and V v = qA u for each charged particle (charge q) in the relation 
Jj! M = — r K^ A A v . The result is a sum over all particles and spin states of the terms 
g 2 7r M!/ . For the electron plasma in particular (2 spin states) 

vr^ = 2e 2 vrH/ e + / g ] . (54) 

It can be extended to the neutrino-plasma system simply by taking the vector po- 
tentials indicated by Eqs. fl37| p8|) namely, —eN 1 + \j2G-pdy for electrons and 
\^2Gpc' v for neutrinos. In this way, one obtains the same relations as (|33| , p4|) 
but with classic theory susceptibilities given by 

< = ^ u [fu + M (55) 

and the tensor above. The differences between this and the field theory results 
( |l^ , |i~3l) only appear, vacuum corrections apart, in the particle propagators and 
parity-violating terms. However, in the limit of frequencies and wavenumbers much 
lower than the electron and neutrino energies, field theory delivers the same results 
as classical kinetic theory. The dispersion relation of the longitudinal waves will be 
analyzed next section . 

4 Waves and plasma instabilities. 

Bingham et al. [|l], |J conceived a mechanism of energy transfer from a neutrino beam 
to an electron plasma in which certain plasma wave modes acquire large growth rates 
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as a result of a neutrino resonant effect. In order that those resonant waves be not 
electron Landau damped the electron plasma has to be in a non-relativistic regime. 
Let ujpiik) designate the frequency of the longitudinal waves as a function of k in a 
plasma without neutrinos. In the case of a Maxwellian distribution uj p i is given by 



wj(fc) = o£ + (56) 

for wavenumbers much smaller than kp>, the Debye wavenumber. uj p is the plasma 
frequency and u 2 = 4Tran e /m e , k^ = 4Tran e /T e for an electron density and tem- 
perature equal to n e and T e , respectively. More important to what follows is that 
for k <C kn (either degenerate or non-degenerate gas) the frequency uj pi does not 
vary much with k and Ttp° M /k 2 is approximately equal to —uj 2 Juj 2 , as long as ui/k is 
much larger than v e x = JT e /m e , the electron thermal velocity. Using this relation 



in Eqs. fljjj , |21| , |23|) one obtains the dispersion relation in the presence of neutrinos 

as 



/2 2 

2 2 (T\ ( - T F C V / ^ /f2 2\2 00 / ctN 

w -Wp,(A;) = ) vtjv , (57) 

up to terms of Gp order (recall that e 2 = Ana), with 

(2# (/ ^ + /p) Ev {k-p v )*-{h*/2y • (58) 

The expressions above put the real impact of the neutrino flux in perspective. 
Keeping only the main factors, 

u -u pl oc w P z — ^ — G F (59) 

clearly indicates that the neutrino contribution is severely suppressed by Gp. The 
only potential exception is a strong neutrino resonance effect. The poles in the 
neutrino propagators represent kinematic conditions for a massless neutrino with 
momentum p to emit or absorb a plasmon with momentum k: (p ± k) 2 = = p 2 . A 
necessary condition for such a Cerenkov process is that k be a space-like vector i.e., 
oj < \k\. That is quite possible for modes that are not Landau damped, \k\ < kn, 
because the Debye wavenumber kp, of a non-relativistic plasma is much larger than 
the plasma frequency to p . 

The frequency and wavenumbers of interest are much smaller than the electron 
and neutrino single particle energies. Therefore, it makes sense to neglect k 2 in front 
of k-p in the neutrino propagators so that 
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oo f <PPv fu + fpk 2 - (k-v u ) 2 



(2tt) 3 (w-jfc.iTJ 



Y2 



This is nothing but the classic theory result contained in the Eqs. ( |53| , |55|), which 
gives also, taking Eq. in account, 

y (2tt) 3 C(J _ fc.^ \ dp dp y 

Together with Eq. (|57D it constitutes the dispersion relation predicted with classi- 
cal kinetic theory. Our result differs from others @, [T^] simply because the forces 
assumed there to account for the weak interactions are different from the Lorentz 
kind of force we derived in section 3. Bingham et al. in particular Jl], U, only con- 
sidered the terms proportional to the electron and neutrino density gradients and 
so obtained a factor of (1 — u 2 /k 2 ) 2 less in the dispersion relation Q. This factor 
comes from the time derivatives of the currents J u and J e in the forces fl39|) and (f4T|): 
J = ukJ°/k 2 for longitudinal waves and the same type of relation holds for the po- 
larization tensor, Eq. (|3]), as a result of plasma isotropy and current conservation. 
Altogether, it makes a factor of 1 — uj 2 /k 2 , one for neutrinos and one for electrons, 
as can be learned from Eqs. ( p3| p4[ ). 



The kinetic instability is analogous to the electron Landau damping |TB[ and 
takes place if the neutrino spectrum crosses from one side to the other of the reso- 
nance (k-v u = u). Then, the integral in Eq. separates into a principal part and 
an imaginary quantity that is evaluated with the so-called Landau prescription 



in this case, with uo — k-v v + i0 + in the denominator. The neutrino contribution to 
the damping rate comes then as 

Glc'v w lT2 2 , 2 f d 3 p u t _ (df v df p \ t 

Hardy and Melrose obtained this result starting from the decay rate of single 
neutrinos into longitudinal photons. However, by the very nature of such calculation 
the full neutrino contribution to the dispersion relation was not derived, which might 
be important to investigate possible reactive instabilities. Before going to that we 
just note that contrary to the interpretation in 0, the factor of (1 — uj 2 /k 2 ) 2 is not 
due to the chiral nature of weak interactions but rather to their current-current 
structure, plasma isotropy and a particularity of longitudinal waves. 

7 W goes as G| and is exceedingly small when compared with the plasma colli- 
sional damping for instance. If however all neutrinos were on the top of the resonance 
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they could generate an hydro dynamic instability. The claim was [I], |J] that this oc- 
curs in the conditions of Supernova neutrino emission causing much larger growth 
rates, proportional to G 2 / 3 rather than G F , or to Gp if electron- ion collisions are 
taken into consideration. The question we raise is, in the end one has to check 
whether or not the entire neutrino flux lies in the resonance i.e., whether \u — k-v u \ 
is confined to the calculated resonance width |7|. If that is not so one falls in the 
Landau damping case and result fl62l) . 

The hydrodynamic limit is obtained by assuming that uj — k-v v is approximately 
constant over the neutrino spectrum in Eq. fl60|) and then solving the dispersion 
relation for a complex u. A solution with positive imaginary part represents a 
reactive instability. From Eqs. (f>7[ |60| ) one writes 

2 2 /A I n v\ ljj2 2\2 k 2 — ik-Vy) 2 . . 

2-Ka\K/ k 2 (u-k-v v ) 2 
where n u means the joint v e and v e flux and we have made c' v = 1. The largest 
growth rates (Im{co>} = —7 > 0) are obtained for resonant modes [k-v v « u p i) with 
magnitudes around 



G F n u G F n e ) 1/3 k 



to: 



But this assumes that • v v covers an interval of width Ak-v u not greater than T. 
The neutrino velocities only spread in direction. Far away from the neutrinosphere 
of radius R v they essentially move in the radial direction yet, the velocity cone at a 
radius r has a finite aperture, 6 V ~ 2R u /r. Hence, 

Ak-v u ^6 u k sin# w ^u;^ |tan^| (65) 

for the resonant modes, where # is the angle k makes with the radial direction. V 
increases with tan 4//3 # so the best chance of satisfying the requirement Y > Ak-v u 
is to have an angle 9 as close as possible to tt/2. It must however not get any 
closer than 7r/2±8 u otherwise Ak-v u becomes larger than u p i and neutrino Landau 
damping cannot be avoided. That puts an upper limit on the ratio V / Ak-v u which 
implies the necessary condition 

1/3 



7i 



f G F n u G F n e } 1 4/3 
I E v m e J " 



Another upper limit and necessary condition comes from that k must be smaller 
than the Debye wavenumber. The necessary condition is 

£ _ f Gf n u G F n e \ 1/3 k D ^ ^ 
I E v m e J 9 U 
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Knowing how small the energies Gpn^ and Gpn e are those conditions look quite 
disfavored by data. Take a neutrino luminosity [fjl| L u = 10 53 ergs/s, E u = lOMeV, 
neutrinosphere radius R u = 30 km and a generous electron density n e = 10 30 cm~ 3 , 
barely compatible with the non-relativistic regime. Recalling that 9 V = 2R u /r and 
n u = L u /47ir 2 one gets 



^ = l - 26 {w^) ■ (68) 

It means that for an electron density as high as 10 30 cm -3 a strong v resonance 
effect could only take place at a radius of 10 14 km or larger! That is clearly absurd. 
The conclusion we draw is the neutrinos are never collimated enough to cause an 
hydrodynamic instability. The resonance width is far too small for that. It slices the 
neutrino velocity cone in two pieces, one above and one below the resonance hv v = u. 
The neutrinos only induce Landau damping with a rate given by the expression 
This was calculated by Hardy and Melrose |1 for an electron density of 10 30 cm -3 
and they found that the growth rates are too small for the time duration of neutrino 
emission in Supernovae. When one tries to push the electron density to increase 
the growth rates one approaches the relativistic regime but then electron Landau 
damping takes over. In addition, the phase space of plasma waves with phase velocity 
less than one (uj < |fc|, necessary for Cerenkov neutrino emission) becomes vanishing 
small thus suppressing the relative factor (1 — uj 2 /k 2 ) 2 in the dispersion relation. On 
the other hand, the collisional damping in a non-relativistic plasma is many orders 
of magnitude larger than the neutrino Landau damping. To conclude, it looks that 
the neutrinos are not capable of transferring any significant energy to the medium 
through plasma wave instabilities. 



5 Conclusions and discussion 

We applied the techniques of finite temperature field theory to study the longitudinal 
modes of the electromagnetic waves in a plasma crossed by an intense neutrino 
flux. It is shown that for an isotropic plasma the electron axial- vector couplings and 
polarization effects are suppressed by Gp in the dispersion relation. Only the electron 
vector couplings (weak and electromagnetic) contribute at G| order. In addition, in 
the limit of frequencies and wavenumbers much smaller than the individual electron 
and neutrino energies, the susceptibility tensors derived from field theory are well 
approximated by the results obtained with classic kinetic theory provided that the 
electroweak interactions are described with the appropriate forces. As a result of 
the vectorial nature of the interactions, the weak forces are of the same type as 
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the Lorentz electromagnetic force. They diverge however from the forces employed 
by other authors Jl|, H @j> XL lYi particular the ponderomotive force of Bingham 
et al. |], |. That expl ains the difference between the dispersion relation derived by 
us and the ones obtained in |4| and |17| using classic kinetic theory. 

In the early papers of Bingham et al. [|l| the neutrinos were treated with a wave 
function obeying a sort of Klein-Gordon equation with an external potential ac- 
counting for the interaction with the medium. The collective effects were attributed 
there to a puzzling phase coherence between the neutrino wave functions. There is 
nothing wrong with that: it simply means that under an external plasma wave of 
wavenumber the neutrino wave function fluctuations share a phase factor, com- 
mon to all neutrinos (e~ lk ' x ), on top of the arbitrary initial phases of each neutrino 
wave function. The problem was rather that in this Klein-Gordon sort of equation 
(which could in principle be derived by squaring the Dirac equation) all the terms 
that might depend on the external potential derivatives or electron 3- vector current 
were completely discarded. That is fine to study problems like neutrino oscillations 
but not for neutrino effects on plasma waves because their very nature concerns the 
variations of the medium densities on the time and length scales of the wave period 
and wavelength. 

The other problem concerns the claimed [p], |[] non-linear effects and plasma 
instabilities induced by resonant neutrinos (k-v u = uo) on a non-relativistic plasma. 
We compared the calculated damping (growth) rate of the plasma waves with the 
phase space occupied by the neutrino flux and concluded that the angular dispersion 
of the neutrino velocity due to the finite size of the neutrinosphere exceeds by far 
the assumed resonance width, at least a factor of 10 8 even in the case of the most 
optimistic growth rate, no matter how large is the distance to the neutrinosphere. 
This is of course due to the severe weakness of the neutrino interactions. As a 
consequence no hydrodynamic (reactive) instabilities can be induced by neutrinos. 
What is left is just neutrino Landau damping, a linear effect proportional to Gp. The 
corresponding growth rates are however too small for the time duration of neutrino 
emission in Supernovae ||. In addition they are overcame by either collisional 
damping in a non-relativistic plasma or by electron Landau damping in a relativistic 
one. In conclusion, wave instabilities induced by neutrinos do not seem to be a viable 
mechanism of substantial energy transfer from neutrinos to a Supernova plasma. 
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